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Abstract. We present metal abundances in 15 damped Lya systems (DLAs) from the Complete Optical and 
Radio Absorption Line System (CORALS) survey, designed to be free from any biasing effects due to extinction 
of QSOs by dust in intervening absorbers. It has long been suggested that such biasing may explain differences 
in metallicity between damped Lya systems and coeval luminous galaxies, and between model predictions and 
observations. We use our measured zinc and chromium abundances (combined with those for five more CORALS 
DLAs from the literature, giving us a very nearly complete sample) to test whether the metallicity and degree 
of dust depletion in CORALS DLAs are significantly different from those of existing, larger, samples of DLAs 
iy-^ ' drawn from magnitude limited, optical surveys. We find that the column density weighted metallicity of CORALS 

DLAs, [((Zn/H) DLA )] = —0.88 ± 0.21 in the redshift interval 1.86 < z a b s < 3.45, is only marginally higher than 
that of a control sample from the recent compilation by Kulkarni et al., [((Zn/H) DLA )] = —1.09 ± 0.10. With the 
present limited statistics this difference is not highly significant. Furthermore, we find no evidence for increased 
dust depletions in CORALS DLAs — their [Cr/Zn] ratios conform to the known trend of increasing depletion 
(decreasing [Cr/Zn]) with increasing metallicity, and we have encountered no cases where Cr is as depleted as 
in local cold interstellar clouds. These results, when combined with the earlier findings of the CORALS survey 
reported by Ellison et al. in 2001, make it difficult to invoke a dust-induced bias to explain the generally low level 
of chemical evolution exhibited by most DLAs. Rather, they indicate that large scale optical QSO surveys give a 
fair census of the population of high redshift absorbers. 
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1. Introduction hydrogen content of the universe available for star forma- 

. . . . tion up to z > 4 (e.g. Storrie-Lombardi & Wolfe 2000 

Spectroscopic studies ol absorption line systems along „ , , u\i\/\n i , ,-, , , , , , , 

, , r . Feroux et al. 2003 ) and arc thought to be the progenitors 

quasar sighthncs are an important source of information c ,. , , , „ T A , r 

n ° . .... . . of the disk galaxies wc sec today. DLAs have therefore 

regarding the chemical evolution history of neutral gas m , , , r , . , , . P , 

° . ^ /T ^ T J . ° been used as a tool for probing the evolution of galaxies, 

the universe. Damped Lya systems (DLAs, defined to have . ln , , . , , ^ .r, , ,. 

on <.,,.,- , , , especially at nign redsnift wficrc direct measurements are 

NHi> 2 x 10 20 atoms cm' 2 : Wolfe et al. 1155511 . which F ' u t T lf , , rmna, t 

v ' — , . . , . „ ^T. — ., . more difncult. Wolfe et al. 11995) found that the comov- 

make up the high column density end of the distribution . , . , . n , n . „ T . ^ 

. , r . J .... . mg mass density m neutral hydrogen m DLAs, ' 'dlAi is 
of absorption line systems, are a particularly interesting . , ., , , , .„ „ , 
z. , , similar to trie mass density m stars at redsnift z — 0, and 
subset ol absorbers to study, ihey dominate the neutral ,, r jj.i_j.-j.ji i j. ,n u 
J therefore proposed that its redshift evolution, di.DLA/dz, 

* This work is based in part on observations collected at could be taken as a global measure of the rate at which 

the European Southern Observatory, Chile (ESO Nos. 69. A- gas in galaxies is converted into stars. 

0051A & 71.A-0067A) and at the W.M. Keck Observatory, 

which is operated as a scientific partnership among the The metallicity evolution of galaxies has also been 

California Institute of Technology, the University of California probed by chemi cal abundanc e measurements m DLAs 

and the National Aeronautics and Space Administration. The (e.g. Pettini et al. 11997b! Prochaska & Wolfe 119991 



Observatory was made possible by the generous financial sup- 12000(1 . which provide the most detailed information on the 
port of the W.M. Keck Foundation. chemical composition of high-z galaxies. If indeed DLAs 
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are representative of the galaxy population at a given red- 
shift, then one would expect that the evolution of their 
neutral hydrogen fraction and metallicity would track that 
of the general galaxy population as a whole. However, no 
strong evolution of either ^dla or metallicity (Zqla) has 
been seen in studies from z ~ 3.5 down to z ~ 0.5 (e.g. 
Rao & Turnshek 120001 Ryan- Weber et al. 120051 Kulkarni 
et al. 12005) . 

Furthermore, it is now well established that DLAs are 
generally at the low end of the metallicity distribution of 
galaxies at redshifts z = 2 — 3 (see, for example, Figure 
32 of Pettini 1200111 . The typical DLA metallicity at this 
epoch is only Zdla — —1-2, or w 1/15 of solar (Pettini et 
al. 1999. Kulkarni et al. 2005), while near-solar metallici- 
ties are common for luminous galaxies detected directly in 
their rest-frame ultraviolet, optical, and far- infrared light 
(e.g. Pettini et al. 120021 Shapley et al. 12001 de Mello et 
al. 12001 Swinbank et al.|200JJ). 

This difference could have a number of causes. 
Selecting galaxies through H I absorption may preferen- 
tially pick out chemically unevolved systems, either be- 
cause galaxies with generally low rates of star formation 
dominate the cross-section (Mo, Mao & White 1H)98[I . or 
because the H I cross-section is largest during a stage prior 
to the onset of star formation. However, direct imaging of 
DLA galaxies shows a varied population of hosts which 
span a range of luminosities and morphological types 
(Boissier, Peroux, & Pettini 120031 Rao et al. 120031 Chen 
& Lanzetta 120031 Weather ley et al. 120051 — see also Ryan- 
Weber et al. I2003[) . Abundance gradients may contribute 
to the difference in metallicity between the outer regions — 
which offer the larger cross-section for absorption — and 
the inner regions of galaxies where star formation activity 
is more prominent (Pettini et al. H994al Chen, Kennicutt, 
& Rauch l200B1 Christe nsen et al. l200Kl Ellison, Kewley, & 
Mallen-Ornclas 2005b). Quantitatively, however, the mag- 
nitude of such gradients has recently been questioned in 
nearby spirals (Bresolin, Garnett, & Kennicutt 2004), and 
remains unknown at high redshifts. 

A third possibility is that damped Lya systems which 
are both metal- and gas-rich do exist, but are system- 
atically underrepresented in current samples drawn from 
magnitude-limited QSO surveys. The hypothesis is that 
even moderate amounts of dust associated with interven- 
ing galaxies may be sufficient to preferentially exclude red- 
dened QSOs from optical surveys, and that the statistics 
of DLAs would accordingly be skewed by such bias against 
dusty absorbers. It is this third possibility which we ad- 
dress in the present paper. 

The idea of dust obscuration of QSOs has a long his- 
tory (e.g. Heisler & Ostriker 1988]) and received observa- 
tional support by the work of Pei, Fall, & Bechtold ( 1991 ) 
who found the spectra of QSOs with DLAs in the redshift 
range 1.77 < z a b s < 2.80 to have statistically steeper con- 
tinuum slopes than those of a control sample. On the basis 
of these results, Fall & Pei l|199 3) proposed that between 
10% and 70% of bright QSOs at z = 3 may have been miss- 
ing from optical samples. This claim is now tempered by 



the recent re- analysis by Murphy & Liske l|2004fl based on 
the much larger compilation of QSO spectra made avail- 
able by the Sloan Digital Sky Survey (SDSS — Stoughton 
et al. l2002J) . From the comparison of 70 QSOs lying behind 
DLAs at 2.0 < z a b s < 4.0) with a control sample which is 
one order of magnitude larger, Murphy & Liske concluded 
that the difference in the continuum slopes a between the 
two sets of spectra is only Aq = —0.04 ± 0.05 (using the 
usual definition of the spectral index, whereby the QSO 
continuum flux is a power law of the form f v oc v a \ corre- 
sponding to a limit on the colour excess due to SMC-like 
dust-reddening of E(B-V) < 0.02 mag (3a). This value 
is significantly lower than Aa = — 0.38± 0.13 reported by 
Pei et al. (BSD). Similarly, Ellison, Hall & Lira l|2()()5a|l 
find E(B — V) < 0.05 mag (3er) from a study of the optical 
to infrared colours of a subsample of CORALS QSOs with 
DLAs in the range 1.8 < z a b s < 3.5. 

However, the SDSS results still do not preclude the 
possibility that highly obscured QSOs may by missing, or 
underrepresented in optical samples (since SDSS QSOs 
are optically selected). Theoretically, such selection ef- 
fects have been appealed to in order to reconcile the pre- 
dictions of hydrodynamic simulations (Cen et al. 2003 
Churches, Nelson & Edmunds 12001 Na gamine, Springel, 
& Hernquist 2004) and galactic chemical evolution mod- 
els (Prantzos & Boissier I2()()()|l with the observations. 
Observationally, their importance is suggested by the ap- 
parent anti-correlation between neutral hydrogen column 
density iV(H i) and metallicity ^dla first pointed out by 
Boisse et al. (fT09"%| . 

It is to assess quantitatively the importance of dust- 
induced bias for the statistics of DLAs that the Complete 
Optical and Radio Absorption Line System (CORALS) 
survey was originally conceived. As the name implies, this 
programme aims at measuring the properties of DLAs in 
a complete sample of QSOs selected at radio wavelengths, 
where dust obscuration is not expected to be an issue. In 
the first stage of the project, Ellison et al. (2001b) iden- 
tified a sample of 22 DLAs from intermediate dispersion 
spectroscopy of all the QSOs (66) with emission redshift 
z cm > 2.2 in the Parkes quarter- Jansky sample of flat- 
spectrum radio sources (Jackson et al. 2002 Hook et al. 
2003). The optical spectra were of sufficient quality to 
measure A^(Hi) in the 22 DLAs, enabling Ellison et al. 
(2001b) to determine both the number density of DLAs 
per unit redshift, n(z), and the corresponding comoving 
mass density of neutral gas, S^dla- The values found, 
n{z) = 0.31±g;gl and logft DL A = -2.59^;^ at a mean 
redshift <z> — 2.37, are higher than the corresponding 
quantities previously determined from optically selected 
QSO samples, but only marginally so. In particular, the 
CORALS survey did not uncover a population of high col- 
umn density (N(Hi) > 10 21 cm~ 2 ) DLAs in front of faint 
QSOs. Within the limitations imposed by the small size of 
their sample, Ellison et al. I|2001bj) concluded that selec- 
tion effects due to intervening dust may at most account 
for an underestimate by a factor of ~ 2 in f^DLA- 
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The H I results alone, however, do not tell us about the 
metal and dust content of CORALS DLAs and whether 
they are higher, on average, than those of the optically 
selected DLAs which have been studied extensively over 
the last fifteen years. These are the questions which we 
explore in the present work. 

Specifically, we have conducted a follow-up programme 
of high resolution spectroscopy of the CORALS DLAs 
aimed at measuring in particular the abundances of zinc 
and chromium. Meyer, Welty & York |[""S89) and Pettini, 
Boksenberg & Hunstead CL989J I1990f) first drew attention 
to the diagnostic value of these two elements. Both are 
iron-peak elements, whose abundances track that of Fe to 
within ± 0.1—0.2 dex in Galactic stars of metallicities from 
solar to about 1/100 of solar (Chen, Nissen & Zhao 120041 
Cayrel et al. 120041 and references therein). In the interstel- 
lar medium of the Milky Way, on the other hand, Zn is 
one of the few elements which show little affinity for dust 
grains, unlike Cr which is usually highly depleted (Savage 
& Sembach 1996). In combination, therefore, these two 
elements can be used to obtain approximate measures of 
the overall degree of metal enrichment, via the [Zn/H] 
ratio, and the fraction of refractory elements locked up 
in solid form, via the [Cr/Zn] ratio. 1 Both elements have 
absorption lines of their dominant ionisation stages in H I 
regions conveniently located at AA2026, 2062 A (Znii) and 
AA2056,2062,2066A (Cm). All of these factors account 
for the fact that Zn n and Cr n absorption lines from DLAs 
have been the target of many studies since the 1990s, even 
though echelle spectrographs on 8-10 m class telescopes 
now afford a more comprehensive assessment of the overall 
chemical composition of QSO absorbers (e.g. Prochaska et 
al. l2001f) . Accordingly, in this paper we focus on the abun- 
dances of Zn and Cr in CORALS DLAs and compare them 
with the large body of such measurements now available 
for optically selected DLAs. 

The paper is organised as follows. In §2 we describe 
the observations, data reduction process and column den- 
sity measurements, while in §3 we list the abundances 
measured in the DLAs in our sample. We compare the 
CORALS Zn and Cr abundances to previous surveys in 
§4 and discuss our findings, together with our conclusions 
in §5. 

2. Observations and Data Reduction 

The CORALS sample of QSOs found to have DLAs with 
absorption redshifts 1.8 < z a b s < z em consists of a 
total of 18 QSOs (22 DLAs; Ellison et al. 12001 b|) . Prior 
to this study, spectra of sufficiently high resolution and 
signal-to-noise ratio (S/N) for abundance determinations 
had already been obtained (either by us or by others) for 
five of these DLAs and have appeared in the literature. 
We therefore describe here only observations made on the 

1 We use the conventional notation whereby [Zn/H] = 
log (Zn/H)- log (Zn/H) . 



remainder of the CORALS sample, details of which are 
collected in Table 1. 

i?-band magnitudes in column (2) have been taken 
from Table 2 of the recent study of Ellison et al. (2005a) 
except where indicated. Emission redshifts z em in column 
(3) have been reproduced directly from Table 3 of Ellison 
et al. (2001bJ, while the values of DLA absorption red- 
shift, z a bs listed in column (4) were determined from the 
observations presented here and are quoted to the preci- 
sion of our measurements. One QSO from the CORALS 
sample, B1251— 407, is too faint (B = 23.7) for high reso- 
lution spectroscopy and no abundance measurements are 
therefore available for the two DLAs identified by Ellison 
et al. in its spectrum. However, we do not expect this small 
gap in our survey to affect the conclusions of the present 
study for two reasons. First, the two DLAs in B1251— 407 
are at higher redshifts than the rest of the CORALS sam- 
ple (z a bs = 3.533 and 3.752). Second, their neutral hydro- 
gen column densities, iV(Hi) = 4 x 10 20 and 2 x 10 20 , are 
at the lower end of the DLA column density distribution. 
Thus, the metal and dust content of these two DLAs, even 
if they turned out to be significantly different from those 
of the rest of the sample, would have only a minor im- 
pact on integrated quantities such as the column density 
weighted metallicity considered in §4. Specifically, if we 
assume that, unlike most DLAs, these two systems have a 
solar abundance of Zn, the column density weighted mean 
metallicity of the CORALS sample would be increased by 
only 0.1 dex. 

2.1. Data Acquisition 

The observations were made between April 2002 and 
March 2004 during several runs on a range of telescopes 
and instruments. The majority of the data were obtained 
in service mode using the Ultraviolet and Visual Echelle 
Spectrograph (UVES; Dekker et al. 2000) on the European 
Southern Observatory Very Large Telescope (VLT). The 
spectra of B0405-331 and B0913+003 were obtained with 
the Echellette Spectrograph and Imager (ESI; Sheinis et 
al. 12*002:) at the cassegrain focus of the Keck II telescope; 
for B0537— 286 we used the Low- Resolution Imaging 
Spectrometer (LRIS; Oke et al. 1995) at the cassegrain 
focus of Keck I. B1230— 101 was observed with the newly- 
commissioned Magellan Inamori Kyocera Echelle spectro- 
graph (MIKE; Bernstein et al.|°2003) on the 6.5m Magellan 
Baade telescope. 

Instrument settings were chosen for each QSO ob- 
servation so as to provide coverage at the wavelengths 
of the redshifted Znii and Cm absorption lines. We 
achieved spectral resolutions varying between 0.14 and 
0.25 A FWHM for the UVES and MIKE data, and of 
~1.16A and 3.87A FWHM for the ESI and LRIS data 
respectively. These correspond to velocity resolutions of 
~ 7^8^, ~13kms _1 , ~48kms _1 and — 140kms _1 for 
the UVES, MIKE, ESI and LRIS spectra respectively. 
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Table 1. Details of the QSOs observed. 



QSO 


B 


2cm 


•2-abs 


Telescope/ 


Resolution 


Integration 


S/N a 






mag. 






Instrument 


(A) 


Time (s) 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(?) 


(8) 


B0335- 


-122 


Zl.D 


"3 /MO 
o.44z 


o i vnnc 

o.i / yyo 


\TJ T /T TAATPQ 
V 1j1 Z/U V Hjo 


n on 
U.zU 


i a qnn 

id ouu 




B0347- 


211 


on n 


Z.y44 


i.y4 ( 


WIT O/TT^TPO 

VIjI z/ U Vhjb 


n i /i 
U.14 


q nnn 

o yuu 





B0405- 


331 


1 O A 

iy.4 


z.o / U 


C^AOQO 

z.ooyoz 


„i r tt /rpoT 

xvecK ii / H/oi 


1.10 


q f?nn 
o OUU 


1 '-. 

io 


B0432- 


-440 


iy.o 


fi/IG 

z.o4y 


o Qnon^ 
z.oUzUD 


ATT T /T T\/"FTQ 
VIjI z/ U V HjD 


n 1 fi 
U.lo 


1 /i /inn 
14 4UU 


1 

1 z 


B0438- 


-436 


on *7 
zU. / 


O OfiQ 
Z.ODO 


O Q/17Q£i 

z.o4 / oO 


VIjI z/ U Vm 


U.lo 


i /i 7nn 
14 /UU 


o ~. 

ZrJ 


B0537- 


286 


ly.o 


o i in 
o.HU 


G7/1£i 

z.y ( 4d 


IVeCK l/ljltlo 


o.o r 


si Ann 
o OUU 


oy 


B0913+003 


Zl. ( 


o.U/4 


z. /4o4 


TV' „ „1 TT /T-inj 

Keck 11/Jibl 


1.21 


o OUU 


Q 


B0933- 


333 


i n T 

iy. i 


o fine 

z.yuo 


O ftOOO 

z.oozo 


VL1 z/UVrjb 


n 1 o 
U.lo 


1 a a nn 
14 4UU 


on 
2U 


B1055- 


301 


19.4 


2.523 


J- . C/UOtJU 


VLT 9 /TTVTCS 

V J-J J- £U 1 \J V J-JkJ 


0.14 


10 800 


15 


B1228- 


113 


20.9 


3.528 


2.19289 


VLT 2/UVES 


0.15 


18 000 


15 


B1230- 


101 


19.7 


2.394 


1.93136 


Magellan I/MIKE 


0.25 


18 000 


11 


B1354- 


107 


19.0 


3.006 


2.50086 


VLT 2/UVES 


0.17 


14400 


21 










& 2.96682 




0.19 




14 


B1418- 


064 


20.4 


3.689 


3.44828 


VLT 2/UVES 


0.23 


14 400 


11 


B2311- 


-373 


19.0 


2.476 


2.18210 


VLT 2/UVES 


0.15 


10 800 


25 



a Signal to noise ratio in the wavelength region encompassing the Znll and Cm absorption lines. 
These two values are taken from Table 3 of Ellison et al. (|2001h>p> . 



2.2. Data Reduction 

Due to the differing nature of the instruments used, the 
spectra were not identically reduced, although the same 
standard steps were incorporated in each case. The re- 
duction of the UVES data was performed with the ESO 
UVES data reduction pipeline, while the other spectra 
were reduced using standard IRAF 2 routines. The two- 
dimensional images were trimmed and the bias level sub- 
tracted using the over-scan regions. Pixel-to-pixel varia- 
tions were corrected for by dividing through by a normal- 
ized flat-field of high S /N that had been produced by co- 
adding several flat-field exposures. Pixels which had been 
affected by cosmic ray hits were corrected for, where pos- 
sible, and the echelle orders or single long-slit spectrum 
were traced. The one-dimensional spectrum was then ex- 
tracted and the sky spectrum subtracted (in the case of 
MIKE taking into account the slit tilt, which is not per- 
pendicular to each order and varies over the CCD). 

The individual sky-subtracted spectra were then wave- 
length calibrated (using comparison lamp spectra typi- 
cally taken immediately following and/or prior to each 
science frame and extracted in the same way as above, ex- 
cept for sky subtraction) and then corrected to a vacuum 
heliocentric wavelength scale. In the next step, the spectra 
of each object were rebinned to a common linear wave- 
length scale, with a bin size close to their original size and 
co-added, rejecting any remaining cosmic rays or bad pix- 
els. The co-addition was performed without any weight- 
ing and the corresponding error spectra were summed in 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation 



Table 2. Transition wavelengths and /-values used in this 
study. 



Species 


Wavelength 


/ 




(A) 




Znii 


2026.1370 


0.501 


Mgi 


2026.4768 


0.113 


Cm 


2056.2569 


0.103 


Cm 


2062.2361 


0.0759 


Znii 


2062.6604 


0.246 


Cm 


2066.1640 


0.0512 



quadrature. Gaussian fits to several emission lines from 
the reference lamp spectra near to the redshifted Zn n and 
Cr II lines were made to estimate the values of the spectral 
resolution (FWHM) listed in Tabic 1. 

2.3. Spectral Line Fitting 

Metal absorption lines associated with the DLA systems in 
each QSO spectrum were initially identified based on the 
redshifts reported by Ellison et al. ( 2001b). The vacuum- 
heliocentric wavelengths measured from our spectra for 
the strongest component in each absorption system were 
then used to determine a more accurate redshift for the 
DLA system, based on the average of the redshifts cal- 
culated from each (non-saturated) line. The spectra were 
then reduced to a rest-frame wavelength scale and cut into 
sections corresponding to ± 2000kms _1 around each ab- 
sorption line. These sections were then normalised to the 
local continuum by division by a spline fit to portions of 
the spectrum judged to be free of absorption. 
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Fig. 1. Sample portions of spectra of CORALS QSOs, showing the Znii A2026.14 and Mgi A2026.48 lines (left), and 
the Cm AA2056.26, 2062.24, 2066.16 and Znii A2062.66 lines (centre) in four DLAs. In the z &hs = 2.50086 DLA in 
the spectrum of B1354— 107 (third panel), we just detect Cm AA2056.26, 2062.24, while the other lines are below our 
detection limit. In the far right panels we have reproduced examples of stronger transitions from other elements which 
were used to determine the velocity structure of the gas in each DLA. The continuous coloured lines show the profile 
fits to the data produced by VPFIT (see text). 



Ion column densities were deduced from the observed 
absorption lines using the line fitting software package 
VPFIT 3 . VPFIT fits multiple Voigt profiles (convolved 
with the instrument profile) to absorption line compo- 
nents. Initial guesses for the column density, redshift 
and Doppler b parameter for each absorption component 
are used as inputs into the program. VPFIT then com- 
putes the x 2 goodness-of-fit (taking into account the er- 
ror on each pixel) and automatically varies the parame- 
ters, iterating until x 2 nas been minimised and the best 
fit has been found. Components were rejected from the 
fit if their Doppler parameters were lower than approx- 
imately half the spectral resolution. For the rest-frame 
wavelengths and /-values of the transitions, which are in- 
puts to VPFIT, we consulted the recent compilation by 
Morton (2003); for reference, values for the Znn and Cm 
multiplets are reproduced in Table 2. 

The first step in the line fitting procedure was to deter- 
mine the multi-component velocity structure of gas in the 
damped system by fitting lines spanning a wide range of /- 



3 VPFIT is available from |http://www.ast .cam. ac.uk/~rfc/| 
vpfit.html 



values, preferably multiplets of Fen and Sin. The param- 
eters of the fit — that is redshift and Doppler parameter — 
were then kept fixed and applied to the generally weaker 
Zn II and Cr II lines. In cases where only some of the Fe n 
and Si n absorption components are detected in Zn II and 
Cm, the total column densities of Zn + and Cr + were ad- 
justed upwards to allow for the unseen components in the 
same proportion as determined for Fen and Sin. The im- 
plicit assumption is that the relative proportions of ab- 
sorbers in each component are the same for all the ions 
considered. While in principle this may not be true, be- 
cause of possible differences in the degree of dust deple- 
tion of different elements, in practice clement ratios are 
found to be remarkably uniform between the multiple ve- 
locity components within a DLA (e.g. Prochaska 2003). 
The correction for unseen components amounted to less 
than ~ 10% of the total column density in all but two 
cases, B0432— 440 and B0933— 333, where the correction 
for unseen components is as much as ~ 50%. Indeed, this 
correction for undetected components was the main reason 
why we used VPFIT to analyse the Zn II and Cr II absorp- 
tion. In all cases considered here, the Znii and Cm lines 
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Table 3. Zn and Cr abundances in the damped Lya systems observed. 



QSO 


2DLA 


JV(H°) 

I 1U Clll j 


iV(Zn + ) 


N(Zn + )/N(U°) 
(W ) 


7V(Cr+) 


N{Cr + )/N(H°) 


A r (Cr + )/AT(Zn + ) 






1 1U Clll 1 


I 1U Llll ) 


(~\ n _9 ^ 
l iU ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


B0335-122 


3.17995 


6.0 ± 1.5 


< 1.8 


< 3.0 


a 


a 




B0347-211 


1.947 


2.0 ± 0.5 


< 2.4 


< 12 


< 7.4 


< 37 




B0405-331 


2.56932 


4 ± 1 


< 5.5 


< 14 


< 21 


< 54 




B0432-440 


2.30205 


6.0 ± 1.5 


< 1.6 


< 2.7 


13 zb 2 


22 zb 6 


> 8.4 


B0438-436 


2.34736 


6.0 ± 1.5 


5.3 ± 0.3 


8.8 ± 2.3 


12 zb 1.4 


19 zb 5 


2.2 zb 0.3 


B0537-286 


2.9746 


2.0 ± 0.5 


< 3.4 


< 17 


< 20 


< 100 




B0913+003 


2.7434 


5.5 ± 1.4 


< 6.6 


< 12 


< 29 


< 53 




B0933-333 


2.6823 


3.0 ± 0.8 


< 0.98 


< 3.3 


< 8.7 


< 29 




B1055-301 


1.90350 


35 ± 9 


8.2 ± 0.5 


2.3 ± 0.6 


41 zb 2 


12 zb 3 


5.0 =b 0.4 


R1 998 — 1 1 ^ 


9 1 Q9SQ 


4 zb 1 


lu.o zn o.o 


9fi 4- 7 


28 zb 2 


RQ 4- 1 8 
uy zn 10 


9 7 + n ^ 

tj. i zn kj.o 


B1230-101 


1.93136 


3.0 ± 0.8 


8.7 ± 1.0 


29 ± 8 


47 zb 4 


160 zb 40 


5.4 zb 0.8 


B1354-107 


2.50086 


2.5 ± 0.6 


< 0.5 


< 2 


7.4 zb 0.8 


29 =b 8 


> 15 


B1354-107 


2.96682 


6.0 ± 1.5 


< 0.85 


< 1.4 


< 6.0 


< 10 




B1418-064 


3.44828 


2.5 ± 0.6 


< 0.96 


< 3.9 


< 5.1 


< 20 




B2311-373 


2.18210 


3.0 ± 0.8 


< 0.66 


< 2.2 


< 3.4 


< 11 





a The Crn lines were not covered by the UVES observations of this QSO as they fell in the wavelength gap between the two 
red-arm CCDs . 



Table 4. Velocity dispersion parameters of the strongest 
component of detected Zn II and Cr II lines 



QSO zdla b 

(kms" 1 ) 



B0432- 


-440 


2.30205 


12 


B0438- 


-436 


2.34736 


8 


B1055- 


301 


1.90350 


16 


B1228- 


-113 


2.19289 


73 


B1230- 


-101 


1.93136 


47 


B1354- 


-107 


2.50086 


12 



are unsaturated and we would have obtained the same val- 
ues of Af(Zn + ) and 7V(Cr + ) — apart from the 'incomplete- 
ness' correction — had we analysed them with the optical 
depth method fHobbs ri974|l which some prefer to profile 
fitting. 

Out of the 15 DLAs listed in Table 1, we report pos- 
itive detections of both the Znii and Cm multiplets in 
four cases (see Table 3). In two additional cases we detect 
Cr II but not Zn II, while in the remaining nine DLAs we 
place upper limits to both iV(Zn + ) and iV(Cr + ). The up- 
per limits are 3er, deduced from the measured S/N ratios 
(listed in column (8) of Table 1) and the velocity spread of 
the absorption implied by stronger lines of Fe n and Si II, 
as explained above. 

Figure 1 shows examples of the spectra in the regions 
of the Znii and Cm multiplets, including three of the 
Zn II and Cr II detections and one case when only Cr n is 



detected. 4 Also shown in the figure are examples of the 
stronger lines used to determine the velocity structure of 
each DLA, as well as the fits generated by VPFIT. Table 
4 lists values of the Doppler b parameter for the detected 
Zn II and Cr II lines; full details of the velocity struc- 
ture of other absorption lines will be presented in a future 
publication (Akerman 2005). 

3. Abundance determinations 

Table 3 lists the column densities of Znii and Cm, or 
upper limits, measured in this study together with the 1 a 
uncertainties in their values. Since the first ions are the 
dominant ionization stages of Zn and Cr in Hi regions, 
we can readily derive the abundances of these two ele- 
ments by dividing the values of 7V(Zn + ) and N(Cr + ) by 
the neutral hydrogen column densities A r (H°) listed in col- 
umn (3) of Table 3 — the results are given in columns (5) 
and (7) respectively. The values of A^(H°) are reproduced 
directly from the compilation by Ellison et al. (2001b) be- 
cause (a) few of our spectra include the Lya absorption 
line and (b) even in cases where we cover the transition, 
the intermediate resolution spectra of Ellison et al. are 
better suited than ours to the determination of A^H ) by 
fitting the damping wings of the Lya absorption line. We 
have adopted throughout a conservative error of ±25% 
to AT(H ), based on the typical accuracy with which such 
measurements have been reported in the literature. 

4 Spectra of all the CORALS QSOs cov- 
ered in the present study are available from 
ftp://ftp.ast.cam.ac.uk/pub/papers/CORALS/ 
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Table 5. Abundances in the CORALS DLA survey. Measurements are given in logarithmic units relative to the solar 
abundances from Lodders 12 + log(Zn/H) = 4.63, 12 + log(Cr/H) = 5.65, 12 + log(Fe/H) = 7.47, 12 + 

log(Si/H) = 7.54. 



QSO z abs 7V(H°) [Zn/H] [Cr/H] [Fe/H] [Si/H] Ref. 

(10 20 cm- 2 ) 



B0335— 122 


3 17995 


6.0 


< 


— 1.16 






—2.61 


—2.56 




B0347— 211 


1.947 


2.0 


< 


—0.54 


< 


—1.08 


< —1.16 






R 04(15— 331 


2 56932 


4.0 


< 


—0.49 


< 


—0.92 


— 1.74 


— 1.40 




B0432— 440 


2 30205 


6.0 


< 


— 1.21 




-1.30 


— 1.45 


— 1.12 




B0438— 436 


2 34736 


6.0 




-0.68 




-1.36 


— 1.30 






B0458— 020 


2 0395 


44.7 




-1.15 




-1.50 


— 1.61 


> —2.10 


1 2 


B0528— 250 


2.141 


8.9 




-1.45 




-1.50 


— 1.57 


— 1.27 


3 4 


B0528— 250 


2.811 


12.9 




-0.47 




-1.11 


— 1.11 


—0.64 


3 4 


B0537-286 


2.9746 


2.0 


< 


-0.40 


< 


-0.65 








B0913+003 


2.7434 


5.5 


< 


-0.55 


< 


-0.92 


-1.59 


-1.47 




B0933-333 


2.6823 


3.0 


< 


-1.12 


< 


-1.19 


-1.54 


-1.22 




B1055-301 


1.90350 


35.0 




-1.26 




-1.58 


-1.57 


-1.13 




B1228-113 


2.19289 


4.0 




-0.22 




-0.81 








B1230-101 


1.93136 


3.0 




-0.17 




-0.45 


-0.63 


-0.24 




B1354-107 


2.50086 


2.5 


< 


-1.32 




-1.18 


-1.25 






B1354-107 


2.96682 


6.0 


< 


-1.48 


< 


-1.65 


-1.54 


-1.31 




B1418-064 


3.44828 


2.5 


< 


-1.04 


< 


-1.34 


-1.72 


-1.48 




B2311-373 


2.18210 


3.0 


< 


-1.29 


< 


-1.60 


-1.70 


-1.52 




B2314-409 


1.8573 


7.9 




-1.01 




-1.17 


-1.29 


-1.03 


5 


B2314-409 


1.8745 


1.3 


< 


-1.18 


< 


-1.59 


-1.85 


-1.86 


5 



References: 1 - Prochaska & Wolfe I0322J, 2 - Prochaska et al. ISuTHV 3 - Centurion et al. ISuTEfl . 4 - Lu et al. ( 1996), 5 - Ellison 
& Lopez jjjQOial . 

In Table 5 we list the abundances of Zn, Cr, Fe, and 
Si, expressed in logarithmic units relative to solar, for 
each of the DLAs in the complete CORALS sample; the 
list includes all the new measurements reported here plus 
the five already available from the literature — references 
are given in the last column of Table 5. The only DLAs 
in the compilation by Ellison et al. Ij2001b(l which are 
missing from our sample are the two systems in front of 
the very faint QSO B1251-407, as explained in §2. We 
have adopted throughout the solar abundances proposed 
by Lodders (2003) in her comprehensive reassessment of 
mctcoritic and photosphcric abundances; those from the 
more recent compilation by Asplund, Grevesse, & Sauval 
(2004) differ by no more than 0.03 dex for the elements of 
interest here. All the Zn and Cr measurements collected in 
Table 5 (and indeed those of the larger comparison sam- 
ple discussed in §4) have been reduced to the same set 
of /-values, as listed in Table 2, except for cases where 
the differences in column density would have amounted to 
less than a few percent. 5 It was not possible, however, to 
ensure the same degree of homogeneity for the Fe n and 
Si II measurements. 



5 Most recent work has used the same set of /-values for the 
Zn II and Cr II multiplets, from the laboratory measurements by 
Bergeson & Lawler ( 1993 ) . Small differences between the values 
quoted by different observers result from rounding errors. 



4. Metallicity and dust in the CORALS survey 

The purpose of this study is to test whether element abun- 
dances in CORALS DLAs are significantly different from 
those of DLAs drawn from optically selected QSO sam- 
ples, and therefore to determine the extent, if any, to which 
previous surveys have been biased against metal-rich, high 
column density absorbers. We consider this question from 
two points of view, by examining first the metallicity dis- 
tributions indicated by the [Zn/H] determinations, and 
then the degree of depletion of refractory elements implied 
by the [Cr/Zn] measures. 

4.1. Comparison of Zn Abundances 

The most extensive, homogeneous, compilation of [Zn/H] 
measurements in DLAs has recently been assembled by 
Kulkarni et al. (2005); it includes 51 detections and 36 
upper limits over the redshift interval z a b s = 0.09 to 3.90. 
(Three of the Znll detections in the Kulkarni et al. data 
set are in common with the CORALS survey and have 
thus not been included in the control sample in the follow- 
ing analysis.) Although larger compilations of DLA abun- 
dance measurements have been published (e.g. Prochaska 
et al. 2003), they bring together data for other elements, 
mostly Si and Fe, in addition to Zn. The improved statis- 
tics afforded by such compilations are offset, in our view, 
by systematic uncertainties due to differing degrees of dust 
depletion and possible nucleosynthctic departures from so- 
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Redshift 

Fig. 2. Comparison of the abundance of Zn measured in DLAs in the CORALS sample (black) and in the control 
sample by Kulkarni et al. 120051 (coloured). Upper limits, corresponding to non-detections of the Znll doublet, are 
indicated by open triangles. The abundance of Zn is plotted on a logarithmic scale relative to the solar value shown 
by the broken line at [Zn/H] = 0.0. 



-1 



-2 - 



Redshift 

Fig. 3. Column density- weighted metallicities of DLAs in 
the CORALS (black) and Kulkarni et al. (coloured) 
samples. For each sample we show two values of [(Zn/H)]; 
the solid squares correspond to the values obtained by 
considering only Znll detections, while the open squares 
(and dashed error bars) are obtained if the upper limits 
to [Zn/H] in Figure 2 are included as if they were detec- 
tions. The squares are plotted at the median redshift of 
the DLAs in each sample. 



lar relative abundances of different elements. Thus, we 
consider it more appropriate to restrict the present anal- 
ysis to the comparison of Zn abundances alone. 

The CORALS [Zn/H] abundances are compared with 
those from the compilation by Kulkarni et al. (2005) in 
Figure 2 (after rescaling the latter to the solar value from 
Lodders 2003). From a visual inspection of the plot we 
conclude that: (a) both the CORALS and the compari- 
son sample of optically selected DLAs are generally metal- 
poor, with typical values of [Zn/H] well below solar; and 
(b) there is a hint that, overall, the CORALS DLAs may 
have marginally higher metallicities, although the consid- 



erable number of upper limits complicates the comparison. 
We now address these points quantitatively. 

The quantity which is of interest for 'cosmic' chemical 
evolution models (Pci & Fall 1995) is the column density- 
weighted mctallicity 



[<(Zn/H) DLA >] = log((Zn/H) DLA ) - log(Zn/H) ( 



where 



((Zn/H) 



DLA/ 



E 

i=l 

E ^(H ), 

1=1 



(1) 



(2) 



which is a measure of the degree of metal enrichment of 
the DLA population as a whole. The summation in eq. 
(2) is over the n DLAs considered in a given sample. The 
present sample is too small to consider subsets of the data 
within different redshift intervals, as done by Pettini et 
al. I|1999|> . Prochaska et al. 1|2003|1 . and Kulkarni et al. 
( 2005). Instead we compute eq. (2) for the CORALS sam- 
ple as a whole, and compare the result with the anal- 
ogous quantity for the Kulkarni et al. (£2)05) sample, 
computed over the same redshift interval spanned by the 
CORALS data, z abs = 1.86 - 3.45 for the full set (57 
DLAs), and z a b s = 1-86 — 2.81 for the Znll detections 
only (27 DLAs). We find [((Zn/H) DLA )] = -0.88 ± 0.21 
for CORALS DLAs, and [((Zn/H) DLA )] = -1.09 ± 0.10 
for the control sample of Kulkarni et al. (2005J. These 
values, which differ at only the 1 a level, were calculated 
considering only the Znll detections. Repeating the calcu- 
lations, but now including the upper limits as if they were 
detections, we obtain [((Zn/H) DLA )] = — 0.87 ± 0.13 and 
— 1.17 ± 0.07 respectively. The errors were estimated us- 
ing bootstrap techniques (Efron & Tibshirani 11993(1 . For 
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Fig. 4. Cr/Zn ratio against Zn abundance for CORALS DLAs (black) and DLAs from the compilation by Kulkarni 
et al. : 2005 coloured). Lower limits, corresponding to cases where Cm lines are detected by Znil lines are not, are 
indicated by open triangles. Abundance ratios are plotted in logarithmic units relative to the solar values shown by 
the broken lines at [Cr/Zn] = 0.0 and [Zn/H] = 0.0. 



each sample (CORALS and the control sample), we con- 
structed random datasets by drawing n times from the 
real dataset with replacement, where n is the number of 
[Zn/H] measurements in each sample. The procedure was 
repeated a million times to build up a distribution of val- 
ues of [((Zn/H) DLA )]; we quote the standard deviation of 
this distribution as our estimate of the error in the quan- 
tity [((Zn/H) DLA )]. 

These results are shown graphically in Figure 3. 
Again we see that the two samples are very similar. The 
marginally higher metallicity of the CORALS sample may 
be real. Alternatively it may be an artifact of small num- 
ber statistics or, in the case where we include upper lim- 
its, it may be due to the higher proportion of upper limits 
skewing the results — over one half of the CORALS mea- 
sures of [Zn/H] are upper limits, compared with one third 
for the control sample of Kulkarni et al. (2005). 

In order to clarify this point, we conducted a statistical 
test based on "survival statistics" (which takes account of 
upper limits) with the program ASURV (LaValley, Isobe 
& Feigelson ll992fl , which implements the statistical meth- 
ods of Feigelson & Nelson 1)1985(1 . Using a Peto-Prentice 
test CLatta lTt)81|l . we tested the null hypothesis that the 
two samples are drawn from the same parent population 
and found the hypothesis to be true at the 90% confidence 
level. We conclude that CORALS DLAs do not exhibit 
significantly different metallicities from those of existing, 
larger, samples of DLAs assembled from optically selected 
QSO surveys. 

4.2. Comparison of Cr/Zn ratios 

An indication of the degree of depletion of refractory ele- 
ments onto dust grains may be obtained from the ratio of 
the abundances of chromium to zinc, as explained in §1. In 



Figure 4 we plot this ratio against the metallicity [Zn/H] 
for each of the CORALS DLAs (coloured) from Table 5 
together with analogous measurements (black) from the 
compilations by Khare et al. (|2004(l and Kulkarni et al. 
( 2005), after rcscaling their values to the same solar abun- 
dances used here. 

Figure 4 shows the trend of increasing Cr depletion 
with increasing metallicity which was previously noted by 
Pettini et al. H1997a[) and shown by Prochaska & Wolfe 
(2002) to be a general feature of refractory elements in 
DLAs. In systems with [Zn/H] ;$ —1.5, [Cr/Zn] is approx- 
imately solar — indicating that there is little dust depletion 
at such low metallicities — while when [Zn/H] > 1, up to 
~ 90% of the Cr can be 'missing' from the gas phase and 
presumably be in solid form. Even so, in none of the DLAs 
do we see the extreme depletions of Cr, by two orders 
of magnitude, commonly measured in cold clouds of the 
Milky Way disk (Savage & Sembach 1996). At the typical 
DLA metallicity, [(Zn/H) DLA ] « —1, approximately 1/2 
to 2/3 of the Cr is in the dust ([Cr/Zn] w -0.3 to -0.5) 
although there is considerable dispersion in the depleted 
fraction Jo- 

All of these facets of the depletion of refractory ele- 
ments in DLAs have recently been discussed by Vladilo 
(2004) who linked them to the (generally early) chemi- 
cal evolution of the galaxies where the absorption systems 
originate. The dependence of /o on [Zn/H] may reflect a 
metallicity dependence of the efficiency of dust formation 
in the ejecta of core-collapse supernovae and in the winds 
of late- type giants. On the other hand, the large scatter 
which accompanies the overall trend is presumably an in- 
dication of how the detailed balance between the processes 
of dust formation, accretion, and destruction is affected by 
the local physical conditions in the ISM. 
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20 20.5 21 21.5 22 

log N(H I) 

Fig. 5. Zn abundance vs. Hi column density for the 20 CORALS DLAs (black) and the 57 DLAs from the compilation 
by Kulkarni et al. (2005; coloured) in the redshift range 1.86 < z < 3.45. Symbols have the same meaning as in Figure 
2. The long-dash line indicates the limit [Zn/H] + logiV(Hi) > 21 beyond which Prantzos & Boissier (2000) noted 
that there are no DLAs. This limit seems to apply to CORALS DLAs too, although it could just be a reflection of the 
small size of our current sample (DLAs with log N (Hi) > 21 are rare). Alternatively, there may be reasons other than 
dust-induced bias to explain this apparent 'zone of avoidance'. Possibly, the cross-section for high iV(Hi) absorption 
decreases with increasing metallicity as gas is consumed by star formation. 




As far as the present work is concerned, it is evident 
from Figure 4 that the depletions of Cr in the CORALS 
sample conform to the overall pattern described above. 
The CORALS DLAs are not extreme in their values of 
/cd nor do they exhibit lower values of [Cr/Zn] at a given 
[Zn/H] than DLAs drawn from optically selected QSO 
samples. 6 If there are any DLAs where the depletions of 
refractory elements approach the high values typical of 
cold clouds in the Milky Way disk, we have not found 
them yet. Based on the data in Figure 4, there is no evi- 
dence that CORALS DLAs should redden the spectra of 
background QSOs any more than a typical DLA from an 
optically selected sample. 

5. Discussion 

The work described in this paper concludes a project be- 
gun six years ago to test the extent to which existing 
DLA samples are biased by dust reddening against gas- 
rich galaxies of high metallicity. The first results, reported 
by Ellison et al. (2001), showed that f^DLA had not been 
significantly underestimated. To that conclusion we now 
add the findings that: 

(1) At redshifts 1.9 < z < 3.5, the metallicity of 
CORALS DLAs, as measured by the [Zn/H] ratio, is only 
marginally higher (at a statistically insignificant level of 

6 Although three out of the four CORALS DLAs with 
[Zn/H] < —1.0 have [Cr/Zn] < —0.5 compared to only five 
out of 16 from the control sample, these two fractions are dif- 
ferent at less than the 1 a significance level even when using 
Poisson statistics (which overestimate the significance of any 
differences between samples when the number of measurements 
is so small). 



only 1 a) than that of DLAs drawn from optically selected 
QSOs— we determine [((Zn/H) DLA )] = -0.88 ± 0.21 for 
CORALS DLAs, and [((Zn/H) DLA )] = -1.09±0.10for the 
control sample of Kulkarni et al. (2005) over the same red- 
shift interval. We note that none of the CORALS DLAs lie 
in the dust- 'forbidden' zone of Prantzos & Boissier (2000) 
where [Zn/H] + log N (Hi) > 21 (see Figure 5). 

(2) The dust-to-metals ratio, as measured by the quan- 
tity [Cr/Zn], exhibits no systematic difference between 
the two samples — we unearthed no evidence to show that 
radio-selected QSOs should be more reddened by inter- 
vening systems than optically selected QSOs. 

These results, together with recent reports of low dust 
extinction in large DLA samples drawn from the Sloan 
Digital Sky Survey (e.g. Murphy & Liske 2004), make 
it increasingly difficult to appeal to dust-induced selec- 
tion effects to explain the observed properties of DLAs. 
It seems unlikely, for example, that the true metallicity of 
DLAs may have been underestimated by as much as a fac- 
tor of five, as recently claimed by Vladilo & Peroux (2005). 
Below the redshift limit of our sample, z < 1.8, the situ- 
ation is less clear. The dust fraction may be higher (e.g. 
Vladilo 2004), but the number density of Mgil absorbers 
from an extension of the CORALS radio-selected sample 
(in the range 0.6 < z < 1.7) is in excellent agreement 
with that from optically selected surveys (Ellison 2005), 
suggesting that dust obscuration is not a problem at these 
lower redshifts either. 

One caveat (often mentioned when considering the re- 
sults of the CORALS survey) is the limited size of the 
CORALS DLA sample: are the 20 DLAs in Table 5 rep- 
resentative of the population as a whole, or are we being 
thrown off course by a statistical fluctuation? Clearly, only 
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future observations of a larger sample of radio selected 
(or possibly X-ray selected) QSOs will settle this issue 
and test, for instance, whether the higher metallicity — 
by 0.21 dex — reported here is a real difference, or possibly 
even an underestimate of the true value of [((Zn/H) DLA )] 
in an unbiased sample of DLAs. For the moment, in the 
absence of better statistics, we can use a Monte-Carlo ap- 
proach to address such concerns, particularly if we wish to 
test for differences at a level as high as the factor of five 
claimed by Vladilo & Peroux H2005|) . 

Specifically, we have investigated the likelihood 
of measuring a column density weighted metallicity 
[((Zn/H) DLA )] < -0.88, as found here, by drawing 20 
DLAs at random from a much larger parent sample of 
DLAs with the column density and metallicity distribu- 
tions proposed by Vladilo & Peroux (2005). These au- 
thors approximated the true H I column density distribu- 
tion with a power law of the form /jv(H i) ^ ^(Hi) - ^, 
and the true metallicity distribution with a Schechter 
function fz oc {Z / Z„) a e~ z l z * (motivated by luminosity- 
met allicity relationship of galaxies in the local universe). 
The least extreme set of parameters among those con- 
sidered by Vladilo & Peroux has j3 — 1.6, a — —0.46 
and log(Z*/Z0) = —0.19; with these values, their column 
density weighted metallicity is [((Zn/H) DLA )] = —0.44. 
Our simulations showed that, drawing 20 DLAs at random 
from this parent population, one would find by chance val- 
ues of [((Zn/H) DLA )] as low as —0.88, or lower, only in five 
cases out of a hundred. 

The absence of a detectable dust-related bias in cur- 
rent magnitude- limited samples of optically selected QSOs 
may appear surprising to some. After all, such a bias has 
long been advocated by theorists to improve the match 
of their models to the properties of DLAs (e.g. Cen et al. 
l2H0H Churches et al.EDISI Nagamine et al. l2"0T)4) . It also 
seemed the natural explanation for the empirical lack of 
DLAs with high column densities of metals highlighted 
by Boisse et al. (|1998i) and still present in the larger sam- 
ple assembled by Kulkarni et al. (2005) and indeed in the 
CORALS sample considered here (see Figure 5). The rea- 
son why most optical QSO samples do not underestimate 
significantly quantities such as IIdla and Zdla compared 
with CORALS was clarified by Ellison et al. (2003) an d 
is related to the shape of the QSO luminosity function. 
Optical QSO surveys will not yield significantly skewed 
DLA statistics provided they reach below the break in 
the QSO luminosity function at B ~ 19. Brighter QSO 
samples, on the other hand, may show a bias from either 
dust extinction or lensing (the two effects would of course 
operate in different directions). 

The main conclusion of the CORALS project so far — 
that there are only minor differences, if any, between DLA 
samples drawn from QSOs surveys at radio and at opti- 
cal wavelengths — can only be regarded as 'good news'. Its 
corollary is that the large data samples being made avail- 
able by major projects such as the Two-degree Field and 
the Sloan Digital Sky Survey afford us an unfettered view 
of the absorber population, although surveys which con- 



centrate only on bright QSOs may not. The challenge is 
now to understand, perhaps with more focused theoretical 
efforts, the rightful place of damped Lya systems within 
the diverse population of galaxies known to inhabit the 
high rcdshift universe. 
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